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Abstract. We determine the asymptotics of the two-point correlation function 
for quantum systems with half-integer spin which show chaotic behaviour in the 
classical limit using a method introduced by Bogomolny and Keating [Phys. Rev. Lett. 
77 (1996) 1472-1475]. For time-reversal invariant systems we obtain the leading terms 
of the two-point correlation function of the Gaussian symplectic ensemble. Special 
attention has to be paid to the role of Kramers' degeneracy. 
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Understanding correlations of energy levels of quantum mechanical systems whose 
classical limit exhibits chaotic motion is one of the major topics in quantum chaos. 
The bridge between quantum mechanics and classical mechanics is provided by the 
Gutzwiller trace formula JIJ which relates the quantum mechanical density of states 
d(E) = J2 n &{E ~ E n ) to a sum over periodic orbits of the corresponding classical 
system, 

d(E) ~ d(E) + JL Y, Y, A lk T,e^ E \ h^O, (1) 

7 keZ\{0} 

where d(E) denotes the mean spectral density (which, by Weyl's law, is of order h~~f for 
systems with / degrees of freedom), and the sum extends over all primitive periodic 
orbits 7 and their repetitions, formally including negative ones. S 7 (E) = § pdx 
denotes the classical action, T 7 is the (primitive) period, T 7 = dS 1 (E)/dE, and the 
amplitude A y k involves topological and stability properties. The conjecture of Bohigas, 
Giannoni and Schmit (BGS) || states that for classically chaotic systems, generically, 
the statistics of energy levels can be modeled by the average behaviour of ensembles 
of random matrices. In the case of systems without spin the relevant ensembles are 
the Gaussian orthogonal and the Gaussian unitary ensemble (GOE/GUE) depending 
on whether the system does or does not possess an antiunitary symmetry like time- 
reversal, see, e.g. ||. In the case of time-reversal invariant systems with half-integer 
spin one also has to deal with the Gaussian symplectic ensemble (GSE). 

The main result in understanding eigenvalue correlations in terms of the underlying 
classical dynamics is due to Berry || . He used the so-called diagonal approximation and 
the Hannay-Ozorio de Almeida sum rule ||, see also ||, to determine the asymptotics 
of the spectral form factor, which is the Fourier transform of the two-point correlation 
function R,2{s), see eq. (0) below. This treatment has recently been generalized to the 
case with half-integer spin J7| using an analogue of the Gutzwiller trace formula which 
includes an additional factor due to spin precession l| H . 



In the case of the GOE and the GUE Bogomolny and Keating flOU , see also [XT], [12 
developed a method for the semiclassical evaluation of ^(s) which yields an additional 
term as compared to the diagonal approximation of the form factor. More precisely, 
their method yields the leading non-oscillatory and the leading oscillatory term of R,2(s) 
as s —>■ oo, whereas the diagonal approximation of the form factor corresponds to the 
leading non-oscillatory term. Recently Haake || proposed a method to adapt this result 
to the case of the GSE. But, surprisingly, although he obtained two terms of the large 
s asymptotics of B,2{s), the method failed to reproduce the leading term. The aim of 
this Letter is to present a slightly different approach to systems with half-integer spin 
which correctly yields both the leading non-oscillatory and the leading oscillatory term. 
Note that |TIJ also includes a remark on GSE asymptotics which, however, is not based 
on semiclassics with spin but on a theorem in random matrix theory and, therefore, is 
not to be confused with the problem addressed here. 

The general method of [K| consists of three main steps. Starting from the 
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observation that trace formulae lead to accurate semiclassical quantization conditions 
Bogomolny and Keating propose to base the semiclassical analysis of spectral 
correlations on such an approximate spectrum. In the course of the calculations they, 
secondly, employ the diagonal approximation as introduced in [[5], || . Finally, they make 
use of the assumption that the oscillating part of the integrated spectral density (i.e. 
the contribution of periodic orbits) behaves like a Gaussian random variable. Here we 
will only briefly sketch the necessary changes to the method of Bogomolny and Keating 
in order to take care of the situation with half-integer spin. For the general formalism 
we refer to |nj |ll|, |I|, || . We will also rely heavily on results of [|7| . 



In order to obtain a simple but efficient semiclassical quantization condition, we first 
integrate ([!]) over the energy E which yields a trace formula for the spectral staircase 
function N(E). Taking into account only orbits up to a time T, which below will be 
chosen of the order of Heisenberg time T# = 2irhd(E), one obtains a truncated spectral 
staircase function Nt{E), and the semiclassical eigenvalues E n (T) can be determined 
from the condition fl3|, \U 



N T {E n {T))=n+ 1 -. (2) 

The trace formula ([1]) can easily be integrated if there is a one-to-one correspondence 
between orbits at different energies (i.e. no bifurcations occur when varying E) and 
from successive integration by parts we see that in leading order in h~ it is sufficient to 
integrate the exponential, i.e 

N T (E) ~ N(E) + J2 E i^ 6 ^ 7 ^' *- ' (3) 

7 fcez\{o} 

|fc|T 7 <T 

where the periodic orbit sum will later be denoted by N^ nc (E). At this point it 
is important to take care of Kramers' degeneracy. If the quantum system, with 
Hamiltonian H, has half-integer spin and is invariant under time-reversal, i.e. [H, T] = 
with T = ia y K, where K is the operator of complex conjugation, then each energy 
eigenvalue has at least multiplicity two. One could now attempt to first calculate the 
correlations for the degenerate spectrum and relate the result to the correlations of the 
non-degenerate spectrum, cf. ||. This strategy is successful for the form factor [/J. 
However, since the truncated spectral straicase function Nj-(E) fails to reproduce sharp 
steps of size two, the quantization condition (0) will not yield degenerate eigenvalues but 
two distinct eigenvalues which both have an additional error. Therefore, we instead take 
Kramers' degeneracy into account at this point by imposing the modified quantization 
condition 

N T (E n (T)) = 2n + l (4) 

which produces a semiclassical spectrum {E n (T)} with Kramers' degeneracy already 
removed. Note that this semiclassical spectrum has mean density d(E)/2, and the 
corresponding Heisenberg time is Tjj = irhd(E). Using the Poisson summation formula, 
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the density of states d(E) of the semiclassical spectrum can be written as 

d(E) := J2 S ( E - E n(T)) = \d T {E) J2(-iy J™ Nt(E) , (5) 

n 1/6Z 

where dx{E) = dNx{E)/dT, see (||). Before we can compare spectral correlations with 
results from random matrix theory (RMT) the spectrum has to be unfolded, i.e. the 
eigenenergies are rescaled such that their mean separation is one. To this end consider 
the spectral interval I = I(E, H) := [E — hw, E + fajj] which contains Nj levels. In the 
semiclassical limit this number can be estimated by Nj ~ 2fuvd/2, where d = d(E), 
i.e. as h — > the length of the interval shrinks to zero but the number of eigenvalues 
contained in I(E,H) goes to infinity, cf. 0. Defining the unfolded eigenvalues by 
x n (T) := E n (T)d/2, the density of states Dt(x), x = Ed/2, of the unfolded spectrum 
{x n (T)} reads 

D T {x) := J2 S ( x ~ x n(T)) = \d T {E) . (6) 

n 

The semiclassical two-point correlation function is defined by 

1 rx+hwd/2 
R 2 (s } I) := -— - / D T (x' + - ) D T (x' - - ) dx' - 1 . (7) 

n^dJx-r^d/2 v 2/ V 2/ 

From the BGS-conjecture we expect that in the semiclassical limit R 2 (s,I) converges 
weakly to the random matrix result (_R5f SE (s) in the case considered here), i.e. 

lim f R 2 (s, I) (f>(s) ds = f R 2 SE (s) <p(s) ds (8) 

h ~->° Jr Jr 

for any smooth test function <p e «S(R). We only aim at providing a periodic orbit 
theory for this relation in the combined limit 

s — > oo , d— > oo and s/d — >0 (9) 

which will allow expansions in s/d. Here d — > oo is a consequence of the semiclassical 
limit and Weyl's law. The asymptotics of the GSE-result reads (see, e.g. JL5|) 

vr cos(2 7 r g ) _ l + fsin(2 7 r g ) 
2 K ' 2 2ns (2ns) 2 ' V ; 

Substituting (^) and (|5|) into (|7|) results in 

R 2 (s,I) = l(d T (E'+i) d T (E>-$ Y^(-ir^< vNT{E ' + ^' NT{E '-^) -1, (11) 

where the brackets denote an average over I(E,h), i.e. (...)e' = 2Tw Ie-Hw •••d-E". By 
a stationary phase argument one easily sees that the terms with v ^ v' are of relative 
order 0(l/rf) in the desired limit ([]), i.e. we have 

R 2 (s,I)~J2rv(s,I) (12) 

with 

r„(a, /) := jJd T (E> + f) d T (E' - f) e ^(^'+f)-iM*'-i)) \ ^ (13) 
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The evaluation of r (s,I) is straight forward and will not be shown here. The result 
corresponds to the usual diagonal approximation of the form factor (cf. |K| |3|) and 
therefore in the present situation reads |7j 

which is the leading non-oscillating contribution of R^ SE (s) as s — > oo ([TTj|) . Here 



'«' indicates that (|Tj) is not just an asymptotic relation but we have also used the 
diagonal approximation which is assumend to be valid in the combined limit (0) . Further 
conditions needed to arrive at (|14]) are hyperbolicity of the translational dynamics and 
the mixing property of the skew product of translational and spin dynamics, see [0 for 
details. We remark that the last condition can be weakened to ergodicity; this result 



will be presented elsewhere [|l6j. Introducing an auxiliary variable s', the contributions 



r„(s, /), v ^ 0, can be written as derivatives, 

(15) 



1 d 2 



r„(a, 1) ~ —— 2 j-j^ J™ {s+S,) $,(5, s>) 
(ykv) 2 osos' 

where we have expanded the smooth part N(E) of Nt{E) in powers of s/d. The 
functions Q u (s, s') are then defined by 

*„(*, s') := / e -K C ^'+t)-^(^-^)) \ . (16) 



The next step lies at the heart of the method of [|l(|. Assuming that the exponent 
of (|16D behaves like a Gaussian random variable G(E') with zero mean we can use the 
identity (exp(iG(E'))) E i = exp((—G 2 (E')/2) E i) and subsequently evaluate the exponent 
in diagonal approximation. This assumption is favoured by a well established conjecture 
on global eigenvalue correlations [|17j]. Employing an expansion in s/d the difference in 



the exponent of (jig) reads 

N^ c (E'+^-N^ c (E'- s i)^J2 J2 ^^^(e^^i-e-^^) .(17) 

7 fcsz\{0} 

\k\Tj<T 

For the square of this expression we again make use of the diagonal approximation, 
which was already needed to evaluate r (s, I), i.e. we only keep contributions of orbits 
with like actions, 

2\ 




(N^ C (E' + |) - N^ C (E' - £))■ 

„ ww^H 1 -™^^))) ■ (18) 

|fc|T 7 <T E , 

where the generic multiplicity of periodic orbits, see, e.g., |7[], will be set to g = 2, since 
we are dealing with time-reversal invariant systems. The remaining sum over periodic 
orbits can be evaluated with the sum rule of 0, which, essentially, is the Hannay-Ozorio 



-teg[^r T ) ■ ( 19 ) 
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de Almeida sum rule ||, additionally taking into account the spin contribution. This 
yields 

(N^(E> + |) - N^(E' - i)) 2 ) « A f 1 ~ COS T ( ^ T/) dT> 

/E n Jo 1 

9 i A + s' 

At this point a short remark is in order. The final result for r u (s,I), v ^ 0, will still 
include the cut-off time T. It has been argued [JTU], ^ that T has to be chosen of the 
order of Heisenberg time, i.e. T = Ciihd, where the constant C has to be determined by 
comparing to the asymptotics of the RMT-result. However, the result for C will depend 
sensitively on the second term in the asymptotic expansion of the cosine-integral in 
(|I9|). Terms of the same order could also arise from non- leading contributions to the 
sum rule, which, unfortunately, are unknown. We are thus unable to give the correct 
sub- leading term in the asymptotic expansion flT9|) . Therefore we conclude that from the 
above considerations we can only obtain the s-dependence of r v (s,I), v ^ 0, but must 
refrain from any discussion about the cut-off time T. Putting together (|15D , ([16|) , (|19|) 
and the Gaussian ansatz we first observe that contributions with \v\ > 1 decay rapidly. 
The leading correction to ([14]) is hence given by ri(s, I) + r_i(s, I), which yields a term 
proportional to cos(27rs)/s. This result is consistent with the leading oscillatory term 
of# 2 GSE (*),cf. (0). 



Summarizing, the method of Bogomolny and Keating |T(J has been applied to 
time-reversal invariant systems with half-integer spin. As in the previously studied 
cases without spin it correctly reproduces the leading non-oscillatory term and the 
s-dependence of the leading oscillatory term of the two-point correlation function as 
s — > oo. Although giving further semiclassical evidence towards the BGS-conjecture 
open questions, as, e.g., a consistent determination of the correct cut-off time T, remain. 
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